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Abstract 

According to the so-called “global polytropic model”, hydrostatic equi¬ 
librium for a planetary system leads to the Lane-Emden differential equa¬ 
tion. Solving this equation in the complex plane, we obtain polytropic 
spherical shells defined by succesive roots of the real part Re(0) of the 
Lane-Emden function 9. Such shells provide hosting orbits for the mem¬ 
bers of a planetary system. Within the framework of the global poly¬ 
tropic model, we study the exoplanet systems HD 40307, n Ara, Kepler-26, 
Kepler-62, and Kepler-275. We give emphasis on comparing our results 
with observations, and with orbit predictions derived from the “general¬ 
ized Titius-Bode relation”. 

Keywords: exoplanets; global polytropic model; planets: orbits; quan¬ 
tized orbits; stars: individual (HD 40307, Ara, Kepler-26, Kepler-62, 
Kepler-275) 


1 Introduction 

Gravitational quantization of orbits in systems of planets or satellites has been 
considered in three recent investigations (i, i, i). In the framework of 
classical mechanics, such systems are assumed to be in hydrostatic equilibrium, 
governed by the Lane-Emden differential equation. We solve this equation in 
the complex plane applying the so-called “complex plane strategy” The 
resulting solution is the complex Lane-Emden function 9. Polytropic spherical 
shells defined by succesive roots of the real part Re(0) seem to provide hosting 
orbits for planets or satellites. In fact, there is only one parameter to be adjusted 
for a particular polytropic configuration defined by 9: the polytropic index n of 
the central body. 

Alternative studies regarding quantization of orbits of planets or satellites 
are cited and discussed in [ij and Q. 
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In the present paper, we emphasize on comparing our results with: (i) cor¬ 
responding observations, and (ii) orbit predictions of undetected exoplanets, 
calculated by the so-called “generalized Titius-Bode relation” (i , Secs. 1-3 
and references therein; the method is analyzed in Sec. 4). Table 2 of shows 
orbit predictions for 29 exoplanet systems; in this table, there are both “interpo¬ 
lated predictions” concerning undetected planets with distances less than that 
of the outermost planet, and “extrapolated predictions” regarding undetected 
planet(s) being next to the outermost detected one(s). For the purpose of ex¬ 
amining, as best we can, the issue (ii) above, we focus on 5 exoplanet systems 
with the larger numbers of predictions: HD 40307 (4-1-1 predictions), /i Ara 
(4-1-1 predictions), Kepler-26 (5-1-1 predictions), Kepler-62 (6-1-1 predictions), 
and Kepler-275 (6-1-1 predictions). 


2 Polytropic Models Simulating Host Stars, and 
Polytropic Shells Hosting Exoplanets 

Preliminary concepts for the global polytropic model can be found in [i| (Sec. 3). 
The cormlex-plane strategy and the complex Lane-Emden function 6 are anal¬ 
ysed in [5| (Sec. 3.1). We will use hereafter the definitions and symbols adopted 

in [I|. 

The real part 6 {^) of 9{^) has a first root -l-iwhich is the radius of 

the central body. Next to ^i, there is a second root ^2 = ^2 +*^o with ^2 > fi, a 
third root ^3 = ^3 + i^o with ^3 > ^ 2 , etc. The polytropic spherical shells S 2 , S 3 , 
..., are defined by the pairs (^ 1 , ^ 2 ), (^ 2 , ^ 3 ), ■ • ■, respectively. A shell Sj can 
be considered as a place appropriate for hosting a planet. For this purpose, the 
most proper orbit radius aj € [Ci-i) is that at which |d| takes its maximum 
value inside Sj, 

max|d[S'j]| = \0{aj -l-iCo)|- (1) 

In the case of two or three planets hosted inside the same shell Sj, there are 
two further proper orbits with radii auj and auj, such that auj < oij < aRj, at 
which |0| becomes equal to its average value inside Sj, 

avg|d[5'j]| = |0(aLj +*6)l = l^(aRj +*Co)|- (2) 


So, planets inside Sj can be hosted on orbits with radii auj, aj, and auj. 

In 0 (Sec. 2), we have developed an algorithm for computing the optimum 
polytropic index riopt of a star, which hosts a number of planets. This algorithm, 
called A[n], can be applied to Np members Pi, P 2 , ■ ■ ■, Pnp of a system with 
Np prescribed distances Ai < A 2 < • • • < A^Vp from the central body. 

The host stars of the exoplanet systems HD 40307 (i, 13 , f§l), /i Ara 


[H, [iil)> Kepler-26 (KOI-2^[T2|, 0), Kepler-62 (KOI-701; 


and Kepler-275 (KOI-1198 


; 0 ) 



are Sun-like stars. We expect therefore that 


the optimum values of the polytropic index n for modeling such stars are about 
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n ^ 3 ( 18|, Sec. 6.1 and references therein; [i|, Sec 3 and references therein). 
Accordingly, we apply the algorithm A[n] to an array {ni} with elements 


4 = 2.500 + 0.001 (z-1), 


i = l, 2, 


1001 . 


(3) 


The 1001 com plex IVPs, counted in this equation, are solved by the Fortran 
code DCRKF54 [l^, which is a Runge-Kutta-Fehlberg code of fourth and fifth 
order modified for solving complex IVPs of high complexity in their ODEs, along 
contours prescribed as continuous chains of straight-line segments (for details 
on the usage of DCRKF54, see [l|, Sec. 4; see also Q, Sec. 3). 

We will hereafter quote the real parts of the complex orbit radii, since these 
quantities have the physical interest in this study. For simplicity, we will drop 
overbars denoting these real parts. 


3 Numerical Results and Discussion 

All astrophysical data used for comparing our results with exoplanet observa¬ 
tions are from the site http://exoplanet.eu. 

In the following tables, the first root of 9, coinciding with the radius of 
the host star, is expressed in both “classical polytropic units” (cpu) — in such 
units, the length unit is equal to the polytropic parameter a ([l|], Eq. (3b)) — 
and solar radii Rq. All other orbit radii are expressed in AU. 


3.1 The HD 40307 system 

Results for the HD 40307 system are shown in Table [TJ For this system, the 
minimum sum of absolute percent errors is 


Amin ^^opt(HD 40307) = 2.531; % = 5, Qc = 6, 


gd = 8, ge = 9, gf = 10, gg = 15 ~ 37.9. 


(4) 


Smaller error is that for g’s distance, ~ 2.2%, and larger one is that for f’s 
distance, ~ 16%. For the most massive planet, d, the error is ~ 3.3%. The 
average error for the computed orbit radii of the 6 planets in HD 40307 is 
~ 6.3%. 

Regarding the large error quoted for f’s distance, an interesting conjecture 
is to associate this distance with the mean-density orbit aRio = 0.2274 AU, 
provided that the maximum-density orbit aio is already occupied by another 
planet not yet observed. In view of this conjecture, the error for f’s distance 
drops to ~ 7.9%, the minimum sum of absolute percent errors drops to ~ 30%, 
and the average error is reduced to ~ 5.1%. 

Table 2 of [J] includes 5 predictions of planet orbits for the HD 40307 system. 
To compare these predictions with orbit radii computed by A[n], we present in 
TableOall internal shells, i.e. shells located before the outermost occupied shell 
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No 15, which seem to be unoccupied according to the up-to-now observations. 
Comparison is also made for the external shell No 17, within which an extrap¬ 
olated planet prediction seems to exist according to 0. The sum of absolute 
percent differences between corresponding estimates is ~ 32%, giving an average 
difference ~ 6.3%. 

The larger difference, ~ 12%, is found for the shell No 12. On this issue, it 
seems interesting to apply again the conjecture discussed above. In particular, 
the prediction P 3 = 0.33 AU can be associated with the orbit aLi 2 = 0.3496 AU, 
leading to a smaller difference ~ 5.9%. However, if, eventually, the average- 
density orbit aLi 2 — 0.35 AU hosts a planet, then the maximum-density orbit 
012 — 0.37 AU should already host a planet. 

3.2 The fi Ara system 

As seen in Table [H the optimum case for the ^ Ara system gives minimum sum 
of absolute percent errors 

Amin^nopt(Ai Ara) = 2.554; = 5, qd = 14, % = 18, qe = 31^ ~ 5.39. (5) 

Smaller error is that for e’s distance, ~ 0.06%, which is also the most massive 
planet; while larger error is that for c’s distance, ~ 3.9%. The average error for 
the computed distances of the 4 planets in fi Ara is ~ 1.3%. 

Regarding predictions of planetary orbits, Table 2 of [3| shows 5 such cases 
for the /i Ara system. Table |4] presents shells providing respective hosting or¬ 
bits computed by A[n]. The sum of absolute percent differences of associated 
estimates is ~ 27%, which leads to an average difference ~ 5.4%. 

The larger discrepancy, ~ 15%, appears in the shell No 8 . Alternatively, we 
can conjecture, as above, that the prediction P 2 = 0.29 AU is associated with 
the orbit ors = 0.2714 AU, deriving a smaller difference ~ 6.4%. Accordingly, 
if, eventually, the average-density orbit ors — 0.27 AU hosts a planet, then the 
maximum-density orbit ag — 0.25 AU should already host a planet. 

3.3 The Kepler-26 (KOI-250) system 

For the Kepler-26 system. Table [S] reveals the optimum case 

Amin ^^opt(Kepler-26) = 2.531; q^ = 5, % = 7, qc = 8 , ge = 18.57. 

( 6 ) 

Here, smaller error is that for b’s distance, ~ 2.9%, while larger one is that 
for e’s distance, ~ 7.2%. The average error in the computed distances of the 4 
planets in Kepler-26 is ~ 4.6%. 

Table 2 in Q shows 6 predictions of planet orbits for the Kepler-26 system. 
In Table[ 6 l we present all unoccupied shells according to up-to-now observations, 
located before the outermost occupied shell No II. The external shell No 12, 
for which an extrapolated planet prediction is made in [^, is also included in 
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this table. The sum of absolute percent differences of respective estimates is 
~ 39%, leading to an average difference ~ 6.5%. 

The larger discrepancy, ~ 16%, arises in the shell No 9. This shell hosts the 
orbit aL 9 — 0.16 AU associated with the prediction P 4 = 0.14 AU, and the orbit 
ag ~ 0.17 AU associated with the prediction P 5 = 0.17 AU. The latter almost 
coincides with its counterpart, while the former is in large discrepancy with 
its counterpart (that quoted above). Since this discrepancy is ^ 3 times the 
average difference, we adjudge this issue in favour of the orbit ang ~ 0.16 AU 
computed by the algorithm A[n], and we thus adopt this orbit in the place of 
the prediction P 4 = 0.14 AU obtained by the Titius-Bode relation. 

The second larger discrepancy, ~ 13%, appears in the shell No 12. We can 
again conjecture that the prediction Pg = 0.27 AU is associated with the orbit 
o:li 2 = 0.2881 AU, giving an appreciably smaller difference ~ 6.7%. In view 
of such a conjecture, if, eventually, the orbit anig — 0.29 AU is occupied by a 
planet, then the orbit ai 2 — 0.30 AU should be already occupied by a planet. 


3.4 The Kepler-62 (KOI-701) system 

For the Kepler-62 system. Tabled gives the optimum case 
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Smaller error is that for b’s distance, ~ 0.08%, while larger error is that for c’s 
distance, ~ 14%. The error for the most massive planet, e, is ~ 0.42%. The 
average error in the computed distances of the 5 planets in Kepler-62 is ~ 4.1%. 

Regarding comparisons with the 7 predicted orbits of the Kepler-62 system 
in Table 2 of [^ , Table | 8 ] reveals that the sum of absolute percent differences of 
respective estimates is ~ 44%, which gives an average difference ~ 6.2%. The 
larger discrepancy, ~ 14%, arises for the orbit aig associated with the prediction 
Pg = 0.55 AU. Instead, we can again conjecture that the prediction Pg is associ¬ 
ated with the orbit ania = 0.5870 AU, resulting in a smaller difference ~ 6.7%. 
Accordingly, if, eventually, the orbit q;li 3 — 0.59 AU is occupied by a planet, 
then the orbit aig ~ 0.63 AU should be also occupied by a planet. In fact, 
we have tacitly invoked this conjecture for the prediction P 5 = 0.33 AU, which 
has been associated with the orbit aRio instead of the orbit aig = 0.2992 AU 
(taking the latter in the place of the former, we find a difference ~ 9.3%, a total 
difference cs 53%, and an average difference ~ 7.5%). 

It is worth noting here that the shells No 5 and No 7, hosting the planets b 
and d on their maximum-density orbits as and ar, do also provide the average- 
density orbits aR 5 and aRv, which are associated with the predictions Pi and 
P 2 , respectively. We can consider this issue as the “complementary case” of 
the conjecture discussed above. In particular, if we assume for a while that 
the planets b and d have not been observed yet, then our eventual conclusion 
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that aR 5 corresponds to Pi and corresponds to P 2 , does also lead to the 
conclusion that 05 and a^ should be already occupied by planets. 


3.5 The Kepler-275 (KOI-1198) system 


For the exoplanet system of Kepler-275, Table |9] gives the optimum case 

Amin ^«opt (Kepler-275) = 2.837; qi, = 4, q'c = 5, gd = 6 ^ ~ 10.5. ( 8 ) 

Smaller error is that for b’s distance, ~ 0.42%, while larger one is that for c’s 
distance, ~ 5.4%, which is also the most massive planet. The average error in 
the computed distances of the 3 planets in Kepler-275 is ^ 3.5%. 

Table 2 of gives 7 predictions for the Kepler-275 system. As seen in 
Table (TUI the sum of absolute percent differences of respective estimates is ~ 
58%, giving an average difference ~ 8.3%. We now focus our attention on two 
cases with large discrepancies. The first case has to do with the prediction 
Pi = 0.03 AU and its associated orbit q;r 2 — 0.022 AU. The resulting large 
discrepancy ~ 25% 3 times the average discrepancy) leads us to adjudge this 

issue in favour of the orbit aR 2 . Second, likewise, the large discrepancy ~ 15% 
2 times the average discrepancy) between the prediction Pq = 0.17 AU and 
its associated orbit aR 5 ~ 0.14 AU turns to be in favour of the latter. Thus we 
adopt the orbits aR 2 and qrs computed by the algorithm A[n] in the place of the 
predictions Pi and Pq obtained by the Titius-Bode relation. Since however Q!R 2 
is an average-density orbit, adopting it as a planet’s orbit prediction, presumes 
that the maximum-density orbit 02 should have priority for hosting a planet. 
It is worth noting here that, for the similar case of cirs, the maximum-density 
orbit as is already occupied by the planet c (Table [S]). 

Finally, it is of particular interest the fact that, in order to associate three 
orbits with the predictions P 2 = 0.04 AU, P 3 = 0.05 AU, and P 4 = 0.06, we 
need to employ all three available orbits in the shell No 3, i.e. Q!L 3 — 0.038 AU, 
as ~ 0.051 AU, and ors ~ 0.061 AU, respectively. Accordingly, any further 
prediction in the interval 0.03 AU, ~ 0.07 AU] should be rejected. 

3.6 Further discussion on the HD 40307 and Kepler-26 
systems 

Tables [T] and [5] show that HD 40307 and Kepler-26 both have optimum poly¬ 
tropic index riopt = 2.531. This means that all respective orbit radii are equal 
when expressed in polytropic units. Equivalently, all respective ratios of orbit 
radii are equal; for example. 



(9) 



O'S / hD 40307 



*^8 / Kepler-26 


= 0.3593, 


( 10 ) 
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etc. 


Next, the orbits and ag are occupied in both systems. In addition, the 
orbit aR ,5 is found to be in satisfactory agreement with respective predictions 
of [3] in both systems. 

Furthermore, the orbits aLe, cte, Olq, and ag in the Kepler-26 system are 
found to be in satisfactory agreement with respective predictions of [^; while 
the orbits and ag in the HD 40307 system are occupied. This means in turn 
that the average-density orbits and aLg could be also predictions for the 
HD 40307 system, but they are “missing” from Table 2 of [4|. 

Relevant worth noting cases are: (i) the orbit ay occupied in the Kepler- 
26 system, and being in satisfactory agreement with respective prediction of 
0 for the HD 40307 system- and (ii) the orbit an occupied in the Kepler-26 
system, and being empty in |^; so, the latter orbit could be a prediction for the 
HD 40307 system, but it is “missing” from Table 2 of Q. 

Finallw the following cases seem to be also “missing predictions” from Ta¬ 
ble 2 of [4|: (i) the orbit aio in Kepler-26, while this orbit is occupied in 
HD 40307; (ii) the orbit aig in Kepler-26, while this orbit is quoted as pre¬ 
diction for HD 40307 in [^; (iii) the orbit ais in Kepler-26, while this orbit is 
occupied in HD 40307; and (iv) the orbit aiy in Kepler-26, while this orbit is 
quoted as prediction for HD 40307 in Q. 


4 Conclusions 

We have examined 5 exoplanetary systems with a total of 22 detected exoplanets. 
Their orbit radii computed by A[n] exhibit an average deviation 4% from 
the corresponding observed distances. In the systems examined, there are 50 
unoccupied internal shells; excluding the 26 internal shells of the /i Ara system, 
this total drops to 24. These 24 internal shells, each providing 3 hosting orbits, 
can accomodate a total of 24 x 3 = 72 exoplanets. We have compared 30 of 
our predictions with respective ones given in Table 2 of 0, estimated by the 
generalized Titius-Bode relation, and we have found an average discrepancy 
^ 6%. In 5 cases, we have applied a conjecture asserting that, eventually, 
hosting orbits are average-density orbits instead of respective maximum-density 
orbits, provided that the latter should be already occupied. Finally, in 3 cases 
with large discrepancies (> 15%), we have adopted our predictions in the place 
of the corresponding predictions of . 
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Table 1: The HD 40307 system: central body i.e. the host star HD 40307, 
and poly tropic spherical shells of the planets b, c, d, e, f, g. For successive 
shells Sj and Sj+i, inner radius of Sj+i is the outer radius of Sj. All radii are 
expressed in AU, except for the host’s radius Percent errors %Ej in the 
computed orbit radii aj are given with respect to the corresponding observed 
radii Aj, %Ej = 100 x \{Aj — aj)\/Aj. Parenthesized signed integers following 
numerical values denote powers of 10. 


Host star HD 40307 - Shell No 

1 



^opt 

Cl (cpu) 

Cl {Rq) 

2.531 (+00) 
5.4317(+00) 
7.16 (-01) 





A 

%E 

b - Shell No 

5 



Inner radius, C4 

Outer radius, Cs 

Orbit radius, Ob = as 

3.7372(-02) 

6.4264(-02) 

4.8999(-02) 

4.68(-02) 

4.70(+00) 

c - Shell No 

6 



Outer radius, Ce 

Orbit radius, Uc = ae 

9.4481(-02) 

8.2084(-02) 

7.99(-02) 

2.73(+00) 

d - Shell No 

8 



Inner radius, Ct 

Outer radius, Cs 

Orbit radius, Od = as 

1.2142(-01) 

1.6372(-01) 

1.3639(-01) 

1.321(-01) 

3.25(+00) 

e - Shell No 

9 



Outer radius, Cg 

Orbit radius, = ag 

2.0736(-01) 

2.0636(-01) 

1 .886(-01) 

9.42(+00) 

f - Shell No 

10 



Outer radius, Cio 

Orbit radius, at = aio 

2.5487(-01) 

2.0838(-01) 

2.47(-01) 

1.56(+01) 

g - Shell No 

15 



Inner radius, Ci4 

Outer radius, Ci 5 

Orbit radius, Og = aig 

5.3513(-01) 

6.3922(-01) 

5.8682(-01) 

6 .(-01) 

2 .20(+00) 



Table 2: The HD 40307 system: unoccupied shells, or unoccupied orbits auj, 
aj, aRj within occupied shells. Pj denote predicted orbit radii of planets not 
yet observed, according to Table 2 of [^; external shells, i.e. shells next to the 
outermost occupied shell, are included in this table only when such predictions 
are available. Percent differences %Dj in the computed orbit radii aj are given 
with respect to the corresponding predicted radii Pj, %Dj = 100x\{Pj—aj)\/Pj. 
Other de tails as in Table [T] _ 




P 

%D 

Shell No 

2 



Inner radius, 

Outer radius, ^2 
Orbit radius, 02 

3.33I0(-03) 

l.I002(-02) 

7.0832(-03) 



Shell No 

3 



Outer radius, ^3 
Orbit radius, 03 

2.0234(-02) 

1.4473(-02) 



Shell No 

4 



Outer radius, ^4 
Orbit radius, 04 

3.7372(-02) 

2.5603(-02) 



Shell No 

5 



Outer radius, ^5 
Orbit radius, ors 

6.4264(-02) 

5.7002(-02) 

6 .(- 02 ) 

5.00(+00) 

Shell No 

7 



Inner radius, 

Outer radius, ^7 
Orbit radius, 07 

9.448I(-02) 

1.2I42(-01) 

1.0609(-01) 

l.l(-Ol) 

3.55(+00) 

Shell No 

11 



Inner radius, ^10 
Outer radius, 

Orbit radius, an 

2.5487(-01) 

3.2477(-01) 

2.8608(-01) 



Shell No 

12 



Outer radius, ^12 
Orbit radius, ai 2 

3.9435(-01) 

3.6973(-01) 

3.3(-01) 

1 . 20 (+ 01 ) 

Shell No 

13 



Outer radius, ^13 
Orbit radius, ai 3 

4.4884(-01) 

4.1694(-01) 

4.4(-01) 

5.24(+00) 

Shell No 

14 



Outer radius, ^14 
Orbit radius, ai 4 

5.3513(-01) 

4.8061(-01) 



External Shell No 

17 



Inner radius, 

Outer radius, ^17 
Orbit radius, ai 7 

7.2715(-01) 

8.2144(-01) 

7.3545(-01) 

7.8(-01) 

5.71(+00) 
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Table 3: The ^ Ara system: central body Si, i.e. the host star /r Ara, and 
polytropic spherical shells of the planets c, d, b, e. Other details as in Table [TJ 


Host star /i Ara - Shell No 

I 



^opt 

6 (cpu) 

a {Rq) 

2.554 (+00) 
5.4898(+00) 
1.245 (+00) 





A 

%E 

c - Shell No 

5 



Inner radius, ^4 

Outer radius, ^5 

Orbit radius, ac = 05 

6.6323(-02) 

l.I495(-01) 

8.7360(-02) 

9.094(-02) 

3.94(+00) 

d - Shell No 

14 



Inner radius, ^13 

Outer radius, ^14 

Orbit radius, aa = 044 

8.2970(-01) 

9.7484(-01) 

9.1279(-01) 

9.210(-01) 

8.91(-01) 

b - Shell No 

18 



Inner radius, ^17 

Outer radius, 

Orbit radius, = ois 

1.45I3(+00) 

1.6443(+00) 

1.5075(+00) 

1.5(+00) 

5.03(-01) 

e - Shell No 

31 



Inner radius, ^30 

Outer radius, ^31 

Orbit radius, ag = 031 

5.0053(+00) 

5.4708(+00) 

5.2381(+00) 

5.235(+00) 

5.84(-02) 
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Table 4: The ^ Ara system: unoccupied shells for which predicted orbit radii 
Pj are given for planets not yet observed, according to Table 2 of 0. Other 
details a s in Table _ 




P 

%D 

Shell No 

6 



Inner radius, ^5 
Outer radius, 

Orbit radius, ag 

I.1495(-0I) 

I.6901(-0I) 

I.4882(-0I) 

I.6(-0I) 

6.99(+00) 

Shell No 

8 



Inner radius, ^7 
Outer radius, 

Orbit radius, ag 

2.1725(-0I) 

2.9575(-0I) 

2.4604(-0I) 

2.9(-0I) 

1.52(+01) 

Shell No 

II 



Inner radius, ^10 
Outer radius, 

Orbit radius, an 

4.6303(-0I) 

5.9148(-0I) 

5.2105(-0I) 

5.2(-0I) 

2.02(-0I) 

Shell No 

24 



Inner radius, ^23 
Outer radius, ^24 
Orbit radius, 024 

2.8301(+00) 

3.1098(+00) 

3.0123(+00) 

2.9(+00) 

3.87(+00) 

External Shell No 

40 



Inner radius, ^39 
Outer radius, ^40 
Orbit radius, 049 

8.9572(+00) 

9.4795(+00) 

9.1042(+00) 

9.17(+00) 

7.18(-01) 
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Table 5: The Kepler-26 system: central body Si, i.e. the host star Kepler-26, 
and polytropic spherical shells of the planets d, b, c, e. Other details as in 
T able [T] _ 


Host star Kepler-26 - Shell No 

^opt 

6 (cpu) 

6 (Rq) 

1 

2.531 (-1-00) 
5.4317(-H00) 
5.9 (-01) 





A 

%E 

d - Shell No 

5 



Inner radius, ^4 

3.0795(-02) 



Outer radius, ^5 

5.2955(-02) 



Orbit radius, ad = 05 

4.0376(-02) 

3.9(-02) 

3.53(-t00) 

b - Shell No 

7 



Inner radius, ^5 

7.7855(-02) 



Outer radius, ^7 

1.0005(-01) 



Orbit radius, Ob = a? 

8.7423(-02) 

8.5(-02) 

2.85(-t00) 

c - Shell No 

8 



Outer radius, 

1.3491(-01) 



Orbit radius, Oc = c^s 

1.1239(-01) 

1.07(-01) 

5.03(-t00) 

e - Shell No 

11 



Inner radius, ^10 

2 . 1002 (- 01 ) 



Outer radius, 

2.6762(-01) 



Orbit radius, Oe = oii 

2.3574(-01) 

2 . 2 (- 01 ) 

7.15(-t00) 
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Table 6 : The Kepler-26 system: unoccupied shells, or unoccupied orbits auj, 
aj, aRj within occupied shells. Pj denote predicted orbit radii of planets not 
yet observed, according to Table 2 of 0. Other details as in Table [21 




P 

%D 

Shell No 

2 



Inner radius, 

Outer radius, ^2 
Orbit radius, 02 

2.7448(-03) 

9.0659(-03) 

5.8367(-03) 



Shell No 

3 



Outer radius, ^3 
Orbit radius, as 

1.6674(-02) 

1.1926(-02) 



Shell No 

4 



Outer radius, ^4 
Orbit radius, a 4 

3.07954(-02) 

2.I0978(-02) 



Shell No 

5 



Outer radius, ^5 
Orbit radius, a^s 

5.2955(-02) 

4.6971(-02) 

5.(-02) 

6.06(-l-00) 

Shell No 

6 



Outer radius, ^5 
Orbit radius, ane 
Orbit radius, ae 

7.7855(-02) 

5.9882(-02) 

6.7639(-02) 

6 .(- 02 ) 

7.(-02) 

1.97(-01) 

3.37(-f00) 

Shell No 

9 



Inner radius, 

Outer radius, ^9 
Orbit radius, ang 
Orbit radius, ag 

1.3491(-01) 

1.7087(-01) 

1.6284(-01) 

1.7005(-01) 

1.4(-01) 

1.7(-01) 

1.63(-1-01) 

2.94(-02) 

Shell No 

10 



Outer radius, ^10 
Orbit radius, aio 

2 . 1002 (- 01 ) 

1.7171(-01) 



External Shell No 

12 



Inner radius, 

Outer radius, ^12 
Orbit radius, ai 2 

2.6762(-01) 

3.2496(-01) 

3.0467(-01) 

2.7(-01) 

1.28(-f01) 


13 













Table 7: The Kepler-62 system: central body Si, i.e. the host star Kepler-62, 
and polytropic spherical shells of the planets b, c, d, e, f. Other details as in 
Table m 


Host star Kepler-62 - Shell No 

1 



^opt 

6 (cpu) 

a (Rg) 

2.773 (-H00) 
6.1063(-H00) 
6.3 (-01) 





A 

%E 

b - Shell No 

5 



Inner radius, ^4 

Outer radius, ^5 

Orbit radius, Ob = 05 

4.3481(-02) 

7.6678(-02) 

5.5342(-02) 

5.53(-02) 

7.65(-02) 

c - Shell No 

6 



Outer radius, 

Orbit radius, Oc = as 

1.1442(-01) 

1.0545(-01) 

9.29(-02) 

1.35(-t01) 

d - Shell No 

7 



Outer radius, ^7 

Orbit radius, Od = 0:7 

1.5633(-01) 

1.2283(-01) 

1 . 2 (- 01 ) 

2.36(-t00) 

e - Shell No 

11 



Inner radius, ^10 

Outer radius, 

Orbit radius, Oe = ctii 

3.5889(-01) 

4.4702(-01) 

4.4511(-01) 

4.27(-01) 

4.24(-01) 

f - Shell No 

14 



Inner radius, ^13 

Outer radius, ^14 

Orbit radius, ai = au 

6.7326(-01) 

7.8344(-01) 

7.1679(-01) 

7.18(-01) 

1.69(-01) 
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Table 8 : The Kepler-62 system: unoccupied shells, or unoccupied orbits auj, 
aj, aRj within occupied shells. Pj denote predicted orbit radii of planets not 
yet observed, according to Table 2 of Other details as in Table [21 




P 

%D 

Shell No 

2 



Inner radius, 

Outer radius, ^2 
Orbit radius, 02 

2.9309(-03) 

1.2620(-02) 

6.2399(-03) 



Shell No 

3 



Outer radius, ^3 
Orbit radius, 03 

2.6848(-02) 

2.3207(-02) 



Shell No 

4 



Outer radius, ^4 
Orbit radius, 04 

4.3481(-02) 

3.0100(-02) 



Shell No 

5 



Outer radius, ^5 
Orbit radius, ttRs 

7.6678(-02) 

6.6406(-02) 

7.(-02) 

5.13(-l-00) 

Shell No 

7 



Inner radius, 

Outer radius, ^7 
Orbit radius, aR 7 

1.1442(-01) 

1.5633(-01) 

1.4004(-01) 

1.5(-01) 

6.64(-H00) 

Shell No 

8 



Outer radius, 

Orbit radius, ag 

2.1787(-01) 

I.8604(-01) 

2 .(- 01 ) 

6.98(-H00) 

Shell No 

9 



Outer radius, ^9 
Orbit radius, ag 

2.7712(-01) 

2.5403(-01) 

2 . 6 (- 01 ) 

2.30(-H00) 

Shell No 

10 



Outer radius, ^10 
Orbit radius, aio 
Orbit radius, orio 

3.5889(-01) 

2.9922(-01) 

3.2838(-01) 

3.3(-01) 

4.90(-01) 

Shell No 

12 



Inner radius, 

Outer radius, ^12 
Orbit radius, ai 2 

4.4702(-01) 

5.4706(-01) 

4.4893(-01) 



Shell No 

13 



Outer radius, ^13 
Orbit radius, ai 3 

6.7326(-01) 

6.2825(-01) 

5.5(-01) 

1.42(-1-01) 

External Shell No 

15 



Inner radius, ^14 
Outer radius, ^15 
Orbit radius, aig 

7.8344(-01) 

9.4463(-01) 

8.4677(-01) 

9.2(-01) 

7.96(-l-00) 


15 
















Table 9: The Kepler-275 system: central body Si, i.e. the host star Kepler-275, 
and polytropic spherical shells of the planets b, c, d. Other details as in TableHJ 


Host star Kepler-275 - Shell No 

^opt 

6 (cpu) 

a (Re) 

1 

2.837 (-hOO) 
6.3115(-H00) 
1.38 (-H00) 





A 

%E 

b - Shell No 

4 



Inner radius, a 

7.1357(-02) 



Outer radius, a 

1.1188(-01) 



Orbit radius, ah = a 4 

9.8416(-02) 

9.8(-02) 

4.24(-01) 

c - Shell No 

5 



Outer radius, a 

1.8278(-01) 



Orbit radius, ac = 0:5 

1.2481(-01) 

1.32(-01) 

5.44(-P00) 

d - Shell No 

6 



Outer radius, a 

3.0822(-01) 



Orbit radius, ad = ckq 

2.3432(-01) 

2.24(-01) 

4.61(-H00) 


Table 10: The Kepler-275 system: unoccupied shells, or unoccupied orbits auj, 
aj, aRj within occupied shells. Pj denote predicted orbit radii of planets not 
yet observed, according to Table 2 of Other details as in Table [21 




P 

%D 

Shell No 

Inner radius, a 
Outer radius, a 
Orbit radius, a 2 
Orbit radius, aR 2 

2 

6.4201(-03) 

2.9258(-02) 

1.3834(-02) 

2.2430(-02) 

3.(-02) 

2.52(-p01) 

Shell No 

Outer radius, a 
Orbit radius, 

3 

7.1357(-02) 

3.8197(-02) 

4.(-02) 

4.51(-H00) 

Orbit radius, 03 

5.1166(-02) 

5.(-02) 

2.33(-H00) 

Orbit radius, ors 

6.0524(-02) 

6 .(- 02 ) 

8.73(-01) 

Shell No 

Outer radius, a 
Orbit radius, and 

4 

1.1188(-01) 

8.2954(-02) 

8 .(- 02 ) 

3.69(-h00) 

Shell No 

Outer radius, a 
Orbit radius, ors 

5 

1.8278(-01) 

1.4480(-01) 

1.7(-01) 

1.48(-H01) 

Shell No 

Outer radius, a 
Orbit radius, aRg 

6 

3.0822(-01) 

2.6992(-01) 

2.9(-01) 

6.92(-l-00) 
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